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For certain crystalline systems, most notably the organic compound EtMe3Sb[Pd(dmit)2]2, ex- 
perimental evidence has accumulated of an insulating state with a high density of gapless neutral 
excitations that produce Fermi- liquid- like power laws in thermodynamic quantities and thermal 
transport. This has been taken as evidence of a fractionalized spin liquid state. In this paper, we 
argue that if the experiments are taken at face value, the most promising spin liquid candidates are 
a Z4 spin liquid with a pseudo-Fermi surface and no broken symmetries, or a Z2 spin-liquid with a 
pseudo-Fermi surface and at least one of the following spontaneously broken: (a) time-reversal and 
inversion, (b) translation, or (c) certain point-group symmetries. We present a solvable model on 
the triangular lattice with an (a) type Z2 spin liquid groundstate. 



The notion of a "spin-liquid phase" - a quantum dis- 
ordered insulating phase which is not adiabatically con- 
nected to a band insulator - has captured the imagina- 
tion of theorists for decades [l|,[2[. In recent years, several 
developments have increased interest in this subject 
including a number of exact results for solvable models 
which prove the existence of spin liquids as theoretically 
stable quantum states of matter [5|-l9|, numerical studies 
[iq|-[i2| of more physically realistic models, increasingly 
sophisticated field theoretic analyses [l, [l3|, and, impor- 
tantly, recent experimental results. Specifically, a num- 
ber of quasi-two-dimensional (2D) insulating materials 
have been found to exhibit highly unusual low tempera- 
ture thermodynamic and transport properties which are 
unlike those expected of conventional phases (3, 16|. 

One of the most notable examples is the organic ma- 
terial EtMe3Sb[Pd(dmit)2]2 (which we will refer to as 
"dmit"), which provides a physical realization of a frus- 
trated spin- 1/2 system on a (anisotropic) triangular lat- 
tice. Although dmit has an odd number of electrons per 
unit cell, the charge response (conductivity) is insulating 
and NMR studies indicate that no magnetic ordering oc- 
curs down to the lowest achievable temperatures, which 
are much smaller than the scale of the characteristic ex- 
change coupling, J ^ 250K [isl]. The specific heat, C, 
uniform susceptibility, and thermal conductivity, f<ixx^ 
exhibit T dependencies consistent with Fermi- liquid- like 
power laws, C ^ k%pT, x ^ PMb' '^xx ^ C', with 
an apparent density of states, p ~ 0.1J~^ per unit cell 
[l6l . [l9| . This is suggestive of the existence of a spin- 
liquid with a charge gap, and neutral, fermionic spinon 
excitations with a non-zero density of states at zero en- 
ergy and an estimated mean free path ~ 0.5 /im [l6| . 
Similar observations in a closely related compound, hz- 
(ET)2Cu2(CN)3 (referred to as kCN),[17] led to an early 
proposal for a candidate spin liquid, with a "pseudo 
Fermi surface" of spinous and an emergent U{1) gauge 
field. [20I Eventuall y it was found that tvCN has a gap to 
mobile excitations, [22| but this does not appear to be the 
case in dmit. 



In this paper, we assume that dmit realizes a fraction- 
alized spin liquid state, and we address the problem of 
identifying the most promising candidates to explain the 
experiments. In particular, we assume that the exper- 
imental claims summarized above can be taken at face 
value, and we also consider the possibility that quenched 
disorder plays no fundamental role. In agreement with 
previous discussions [2o|, we argue that under these as- 
sumptions the U{1) spin liquid is not a viable candidate. 

Firstly, we introduce a Z2 spin Hquid phase with 
a stable pseudo-Fermi surface and spontaneously bro- 
ken time-reversal and inversion symmetry. Other at- 
tractive candidates with pseudo-Fermi surfaces which 
are only marginally unstable are a Z4 spin liquid with 
no necessary broken symmetries, or a Z2 spin liquid 
which spontaneously breaks translation symmetry [21|, 
or certain point-group symmetries. In particular, we 
have studied the magnetic field response of these spin 
liquids, and conclude that these states are compati- 
ble with existing magneto-thermal transport measure- 
ments. If weak quenched disorder does play a funda- 
mental role, then other spin liquid states may also be 
viable candidates. (231 [2^ The broken symmetries of the 
Z2 spin liquids with a pseudo- Fermi surface imply at least 
one thermal phase transition. 

''Parent" U{1) spin liquid - A physically motivated 
starting point for the discussion of spin liquids Q be- 
gins with the representation of the spin operator as 
S = ^^r^, where ^ is a two-component spinor field 
(corresponding to the two polarizations of a spin- 1/2 
fermion), and are the Pauli matrices. This leads to a 
minimal model of a sea of spinous coupled to an emergent 
U{1) gauge field, given by the Euclidean Lagrangian: 



dt — iao — e(— iV — a) 
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e = eol+X^a "^a^a is a 2 X 2 matrix; time reversal symme- 
try implies that eo(^) = eo{—k) and Ja{k) = —Ja{—k), 
and, in the absence of spin-orbit coupling, spin-rotational 
symmetry implies Jq, = 0. Here, k signifies the Bloch 
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wave- vector and, where there are multiple spinon bands, 
it is implicitly assumed to include a band index, /^j, is 
the field strength corresponding to the emergent gauge 
field, a^, and . . . represents four-fermion and higher or- 
der interaction terms, subject to the constraint of spinon 
number conservation which is implied by the emergent 
gauge symmetry. L is proposed to describe physics at 
energies small compared to the charge-gap, so the re- 
maining degrees of freedom carry zero electro-magnetic 
charge. 

The weak-coupling {g 0) fixed point is unstable. 
One possible result is a strong-coupling non- Fermi-liquid 
phase which does not break any symmetries, does not 




have any well defined c 
the Fermi surface. |2. 



juasiparticles, but which preserves 
However, assuming this to be a 
stable phase, it cannot be responsible for the experimen- 
tal claims summarized above. It would exhibit power- 
laws (for example, C ^ T^/^) that are substantially dif- 
ferent from those of a Fermi-liquid, unless a broad in- 
termediate finite-temperature regime is assumed that is 
governed by the unstable {g 0) fixed point. Given that 
the instability is strongly relevant and there are no nat- 
urally small parameters in the problem, we find such a 
broad intermediate regime to be unlikely. Additional is- 
sues that cast serious doubt on the U{1) spin liquid below 
1 K on the basis of the thermal Hall data were presented 
inRefs. Eliel. 

Breaking U{1) to Z2 - The gauge symmetry in Eq. (1) 
can be spontaneously broken due to pairing of spinous, 
gapping the gauge fluctuations and leaving a residual Z2 
gauge symmetry. Because a Z2 gauge theory has no finite 
temperature transition in two dimensions, spinon pairing 
defines a crossover rather than a phase transition. The 
discreteness of the residual gauge symmetry implies that 
there are gapped, vortex-like excitations ( "visons" |27| ) 
which carry half a quantum of gauge flux. Visons are 
time-reversal invariant objects - a point to which we will 
return below. Now, the spinon number is only conserved 
mod 2, so new terms in the effective field theory are per- 
mitted. 

At energies low compared to the vison creation en- 
ergy the resulting effective field theory is of the form: 
L = ^t[5^ + i(-^V) + S{-iV)]^ + ^tA(-zV)^t + 

^A^(— iV)^ + • • • . The induced changes in the disper- 
sion, S{k)^ can be absorbed into a redefined spinon "band 
structure," e(A:) = e(^) -\- S{k). 

One possibility is that pairing fully gaps the Fermi sur- 
face, leading to a class of gapped, topologically ordered, 
spin liquid states, which are also ground states of many 
known solvable microscopic models [HI . These states ap- 
pear to be irrelevant for describing the physics of dmit, 
which appears to have delocalized gapless modes. 

Nodal Z2 spin liquids - If the resulting lines of gap 
nodes intersect the spinon Fermi surface, the result is 
a nodal spin liquid: the gauge fluctuations are gapped. 



FIG. 1: (a) Loop order on the triangular lattice where (j) 
represents the phase accumulated in circling the plaquette. 
This leads to the dispersion relation e{k) = — t[cos(A;a;) + 
cos(A;+) + cos(/c-)] — S[sm{kx) + sin(A;+) + sin(/c_)] where 
k± = -ka:/2 ± V3ky/2 and 6/t = tan((/)/3). (b) The Fermi 
surface with (j) — 7t/2 (solid line) and the same curve with 
k ^ —k (dashed line), (c) The dispersion for ky — in the 
presence of a small pairing term. 



but the spinon spectrum consists of a number of gapless 
Dirac-like nodal points. Dirac nodes are known |28| to be 
perturbatively stable as long as time-reversal symmetry 
is preserved. There are various exactly solvable models 
that exhibit such a stable nodal Z2 spin liquid state 0, Hj 
While this class of states has gapless spinon excita- 
tions, it has a vanishing density of states at zero energy, 
and so, it is not a candidate to explain the experiments 
if the effects of disorder are negligible. It has been noted, 
however, that weak disorder broadens the nodes, result- 
ing in a constant density of states at zero energy propor- 
tional to the spinon scattering rate, 1/r. Thus, only if 
the disorder plays a significant role in the thermodynam- 
ics can a nodal spin-liquid be a candidate to explain the 
experiments. 

Z2 spin-liquids with a pseudo-Fermi surface - In the 
U{1) spinon Fermi surface state, arbitrary pairing terms 
cannot appear in the theory as they are forbidden by the 
U{1) gauge symmetry. However once the U{1) is bro- 
ken to Z2, the spinon number is only conserved mod 2, 
so any pairing terms can generically appear. A small 
s-wave pairing term, which is allowed by symmetry, can 
thus open a gap everywhere, while p-wave or d-wave pair- 
ing terms will gap the Fermi surface everywhere but at 
isolated nodal points. Indeed, in Ref. 0, the 'T matrix 
model" was found to have a Z2 spin liquid groundstate 
which, under fine-tuned conditions, has a pseudo-Fermi 
surface. However, any change in the parameters away 
from this multicritical point results in a nodal spin-liquid. 

The instability inherent in a state with a spinon Fermi 
surface with a Z2 gauge field is removed if the system 
in question breaks both time-reversal and inversion sym- 
metry. In this case, the degeneracy of states at k and 
— ^ is lifted. (See Fig. 1.) If we draw the Fermi sur- 
face corresponding to e(— /c), (dashed line in Fig. lb) a 
small pairing term will open gaps only in the neighbor- 
hood of the points at which the two copies of the Fermi 
surface happen to cross. The rest of the Fermi surface is 
perturbatively stable! 

Therefore, if both time-reversal and inversion symme- 
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try are broken, it is possible to stabilize a Z2 spin liquid 
with a pseudo- Fermi surface. Such a state was found 
in [qI in a F-matrix model on the Kagome lattice with 
explicit time-reversal and inversion symmetry breaking. 
A related result was obtained in another context: in 
[28| it was shown that a Fermi surface occurs in a state 
with coexisting d-wave superconducting and orbital loop 
orders, [m An example of such a state on a triangular 
lattice is shown in Fig. 1. A subtle point is the possibil- 
ity that the spinon band structure violates time-reversal 
and inversion, while physic al (g au^e invariant) quantities 
preserve these symmetries [30|. However such states are 
believed to be unstable to translation symmetry breaking 
that gaps the Fermi surface. 

To affirm the possibility of a stable Z2 state with 
pseudo- Fermi surface, we introduce (in the Supplemental 
Material) an exactly solvable spin- 7/ 2 version of the F- 
matrix model on the triangular lattice. The model itself 
is not inversion symmetric. We find that for a range of 
parameters, the ground-state spontaneously breaks time- 
reversal symmetry and supports a stable pseudo Fermi 
surface coupled to a Z2 gauge field. 

Another possibility is a pair density wave (PDW) state 
for the spinous. A specific version of this was proposed 



2l| . The essential feature [3l| of this state which pre- 



vents it from fully gapping the pseudo-Fermi surface is 
that the gap parameter changes sign under translation 
by 1/2 the PDW period, so that the spatial averaged 
gap vanishes. In the special case where the period is two 
lattice spacings, the PDW does not actually break trans- 
lation symmetry, since translation by one lattice constant 
is equivalent to a uniform gauge transformation. On the 
triangular lattice such a state does, however, break the 
Ce rotational symmetry. (However, on a square lattice, 
a period 2 PDW state could occur without breaking any 
symmetries.) In these cases, the pseudo- Fermi surface 
has a marginally relevant Cooper instability. 

Z4 s^in liquids - We can also imagine cases where 
quartets condense but not pairs, leaving a Z4 gauge the- 
ory with a richer collection of vortex-like modes. While 
quartet condensation is relatively unnatural for electrons, 
with their strong repulsive interactions, there is no par- 
ticular reason to rule out condensation of higher multi- 
plets in the case of spinous. It is, for example, possible 
to construct model electron problems with strong, spin- 
dependent attractions which exhibit a charge 4e super- 
conducting phase. |32j] 

A charge 4e condensate can also occur if a pair- 
density wave state [s^ (PDW) is quantum or ther- 
mally melted by the proliferation of double dislocations, 
and similar considerations apply to FFLO states with 
no phase twist. [35| Thermal melting of a PDW can pro- 
ceed in several ways, but if it occurs by the proliferation 
of double dislocations, the result is a spatially uniform 
charge 4e superconductor. While the theory of the quan- 
tum melting has not been worked out, presumably the 



proliferation of double dislocation loops would similarly 
lead to a charge 4e superconducting phase. Moreover, 
by analogy with the similar problem of stripe melting 
considered in [36|, we expect that the Fermi surface of 
gapless quasiparticles will survive this transition intact. 
In the spin-liquid context, such a state would be a Z4 
spin-liquid with a pseudo Fermi surface. 

A state with a Z4 condensate is exotic, in the sense that 
there is no simple (quadratic in fermion operators) mean- 
field description possible. Certainly, there are no gap- 
less gauge fluctuations and the "visons" necessarily break 
time-reversal symmetry - a vison which carries 0o /4 flux 
is distinct from its time-reversed anti- vison which carries 
— (/)o/4. If the quartet binding energy is large enough, 
the spinon spectrum will be fully gapped. However, the 
melted PDW provides a tangible example of a Z4 liq- 
uid with a Fermi surface. Moreover, such a spin liquid 
would be at most marginally unstable {i.e. that there is 
a Cooper instability to a spinon-paired state.) Any po- 
tentially relevant pair-field perturbation is forbidden by 
the residual Z4 gauge symmetry. A Z4 spin liquid is an 
attractive candidate for accounting for the experiments - 
in contrast to the Z2 spin Hquid, it does not necessitate 
time-reversal and inversion symmetry breaking to ensure 
the stability of the spinon- Fermi surface over a broad in- 
termediate energy scale. 

Response to magnetic field - Experiments on dmit have 
reported no measurable thermal Hall angle up to an ap- 
plied field of 12 T [l6|, and an interesting upturn in 
tvxx{T = 0) starting at an applied field of 2 T. 

The orbital coupling to charge fluctuations leads to a 
linear coupling between the magnetic field and the spin 
chiralitv. |37| The spin chirality is proportional to the 
magnetic flux of the emergent gauge field. Q However, 
the visons of Z2 spin hquids typically are even under 
time-reversal, due to tunneling between 0o/2 and — (/)o/2 
vortices. For strong enough tunneling, vortices are there- 
fore not induced by a magnetic field, even if it is larger 
than the vison gap. If the tunneling amplitude is suf- 
ficiently small, an applied magnetic field could mix the 
symmetric and antisymmetric states, stabilizing vortices 
relative to anti- vortices. When the Z2 gauge field is cou- 
pled to gapless spinous, this tunneling may be further 
suppressed, and could even vanish [3 81] . 

If the effective penetration length of the spinon con- 
densate is small (i.e. if it forms a Type I superconduc- 
tor), then we expect that a magnetic field will not in- 
duce a finite density of vortices, and therefore there is no 
mechanism for modifying the thermal Hall effect. If the 
spinon condensate forms a Type II superconductor, then 
a magnetic-field dependent thermal Hall effect is possi- 
ble in principle, although at present we do not have a 
theoretical estimate of the magnitude. 

Assuming the Type I scenario, the only response of 
the spin liquid to the magnetic field is through the Zee- 
man coupling. For the pseudo-Fermi surface states dis- 



4 



cussed in this paper, this wih change the density of 
states and wih not modify the Hah response. How- 
ever, since the spin hquid with a pseudo- Fermi sur- 
face must spontaneously break time-reversal and inver- 
sion symmetry, it should have a zero-field anomalous 
thermal Hall response. In the clean limit, the only 
contribution to the thermal Hall conductivity comes 
fromjsil the Berry curvature of the Bloch states of the 

spinous: K:^y = ^^f^^ / d'^kfxy{k)n{k), where fxy{k) 
is the Berry curvature, and n{k) is the occupation num- 
ber of the partially filled bands. Generically, we ex- 
pect: ^ / (Pkfxy{k)n{k) ~ 1. Therefore, we estimate: 

i^xy/T ^ ^ 10"^^ W/K^ . In order to compare 

with the experimental results, consider Hxy/T per layer, 
where the layer spacing is d ~ 1.5 nm: Kxy/Td ~ 6 • 10~^ 
W/K^m. Using the zero temperature intercept of the 
longitudinal thermal conductivity, limT^o t^xx/Td = 0.2 
W/K^m, we predict a Hall angle id^nOH = i^xy/i^xx ^ 
0.003. The experimental error bars on the Hall angle 
were reported to be ^ 0.05 at 0.23 K, ^ 0.02 at 0.70 K, 
and 0.003 at 1 K and 12 T [l6|. Our prediction is over 
an order of magnitude less than the experimental error 
bars on the low temperature data and therefore appears 
to be consistent with experiment. 

The other observed feature in the magneto-thermal 
transport is a rapid upturn in k,xx{T = 0) at 2 T. One 
possible explanation is that the tunnel splitting in the 
vison state is small, so that the gap to vortex formation 
is approximately li^, and hence magnetic fields in ex- 
cess of 1 T lead to vortex proliferation. (There is also 
a peak in tixx{T)/T at T ^ IK.) However, the vortex 
state would be expected to exhibit a thermal Hall effect 
due to scattering of the spinous from vortices, and has 
not been observed experimentally. One possibility is that 
the experimental error bars are still too large to observe 
this effect. An alternative possibility is that the upturn 
in Kxx/T is simply due to variations of the density of 
states in the spinon band structure. 

Since the Z4 spin liquid with a pseudo-Fermi surface 
need not break time-reversal symmetry, its anomalous 
thermal Hall response can be exactly zero. The structure 
of its vortices is richer, and in particular the Z4 visons are 
not all time-reversal symmetric and will therefore couple 
linearly to a magnetic field, leading to the possibility of 
a non-zero thermal Hall effect. However, the same quan- 
titive issues discussed in the Z2 context apply here, as 
well. 

Discussion - So far, we have framed the discussion in 
terms of the stable spin liquids derived form the "parent" 
U(l) spin liquid. More generally, all known stable spin 
liquids can be understood in terms of a gauge field with 
gauge group G coupled to matter fields. If the spinous 
have a Fermi surface which is at most marginally unsta- 
ble, then presumably the gauge group must be discrete. 

Nodal fermions, with dispersion e{k) ~ /c" are another 



possibility, although for the familiar case of a = 1 the 
effects of disorder must be invoked to account for a non- 
zero density of states. Moreover, in this case the thermal 
conductivity is theoretically expected to be "universal" 
in the sense that it does not depend on r, although it 
does depend on the anisotropy of the nodal dispersion; 
as we will show in the Supplemental material, to account 
for the magnitude of the observed thermal conductivity, 
one must assume an extremely large anisotropy corre- 
sponding to a pairing scale that is roughly 500 times 
smaller than J. In contrast, neither of the proposed 
states with pseudo-Fermi surfaces require the existence 
of such an unnaturally small pairing scale. The case of 
a = 2 ("quadratic band touching") does not seem to 
have been considered previously; this would produce a 
finite density of states, although at the expense of ren- 
dering the state marginally unstable in the presence of 
interactions (ioj . 

As we have discussed, the pseudo- Fermi surface is com- 
pletely stable in the Z2 case only if time-reversal and in- 
version are broken, and marginally unstable if only trans- 
lation and/or rotational symmetry is broken, (isj Other- 
wise, we require a larger gauge group, such as Z4 Q. 
An interesting possibility is given in [24], which consists 
of both a pseudo- Fermi surface of spin-1 fermionic exci- 
tations and a node with dispersion e{k) ~ k^ . This state 
also breaks time-reversal and inversion, although it also 
requires weak disorder to exhibit power-laws consistent 
with experiment. 

One experimental signature of the chiral Z2 states is 
that the breaking of time-reversal symmetry must occur 
at a finite T transition (presumably in the Ising univer- 
sality class) [ill , and such a transition should be observ- 
able in any thermodynamic quantity. Existing specific 
data show no sharp anomaly, but other thermodynamic 
quantities, such as elastic moduli, might be more sensi- 
tive. In the time-reversal breaking phase, various anoma- 
lous response functions should be non-zero. For exam- 
ple, as we have discussed, one expects a non-zero anoma- 
lous thermal Hall effect, although it may be small. One 
also expects a spontaneous Kerr effect i.e. an anoma- 
lous Hall effect, which can be measured with exquisite 
sensitivity. [42^ When spin-orbit coupling is taken into ac- 
count, there should be small magnetic fields which might 
be detectable in NMR or /xSR, although since these fields 
are proportional both to the magnitude of the time- 
reversal symmetry breaking order parameter and to the 
spin-orbit coupling, they may be quite small. 

The marginal instability of the time-reversal invariant 
spin liquids with pseudo-Fermi surfaces has possibly in- 
teresting experimental implications; it could account for 
the existence of small energy scales in the problem which 
could depend sensitively on minor differences between 
materials. For instance, it might account for the low 
temperature gap inferred from thermal transport in a^CN 
(22} . Moreover, the Z4 scenario does not require any fi- 



5 



nite temperature phase transitions. 
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Supplemental Materials 

A. Gamma matrix model with a Z2 spin liquid phase 
which spontaneously breaks TR symmetry and has a 
spinon Fermi surface 

As a proof of principle, we consider an exactly solv- 
able spin-7/2 model on the triangular lattice to illustrate 
a Z2 spin liquid with a stable pseudo-Fermi surface with 
breaking time reversal and inversion symmetries. Follow- 
ing the spirit of the original F-matrix models [8[, we first 
write down the F-matrices in terms of spin-7/2 operators: 



— j^Q^S^ ^ a — 1, • • • ,7, 



(Al) 



where are the spin-7/2 operators and /^^^ are a set of 
numbers chosen so as to satisfy the Clifford algebra of the 
Gamma matrices: {F^, F^} = 2(5^^. It is known that /^^^ 
can always to be chosen symmetric: /^^^ = /^^^ = Z^^^. 
We denote F^'^ = [F^, r^]/(2z). Both F^ and F^'^ are odd 
under time reversal. 

In terms of the F-matrices defined above, we consider 
the following Hamiltonian: 

3 

TT X^XT^ \l r2a-lr2a , 7/ T^2a-l,7T^2a,7 



7 -p2a — l-p2Q!,7 I 7/ -p2Q! — l,7-p2Q! , 7 -p7 
Joc^ i I i+e« ^ ^a^i ^ i+e. + -^71 i 



{PP') 



(A2) 



where Jq,, J^, Jq,, and are coupling parameters, Wp 
are the plaquette operators defined for each triangular 
plaquette, {pp') denotes two edge-sharing triangular pla- 
quettes. The plaquette operator Wp is defined as 



TT/^ _pl6p23 p45 
^^p—^i J- 2+ei+e2 



for up-triangles and 



t^16t^23 t^' 
^ i ^ i — ei ^ 



45 



(A3) 



(A4) 




FIG. 2: (a) The schematic representation of the triangular 
lattice, (b) The Fermi surface of Ck in Eq. (|A8P with Ja — 
{0.3,0.3,0.3}, J4 = {0.3,0.3,0.3}, = {1.0,1.0,1.0}, = 
{-0.8, -1.0, - 1.2} and J7 = 0.5. (c) The Fermi surface of 
^_(k) in Eq. (|AlTl) . 



for down-triangular plaquettes. Note that the Hamil- 
tonian is time reversal symmetric but explicitly breaks 
inversion symmetry of the triangular lattice. 

We represent the Gamma matrix operators in terms of 
bilinear Majorana fermions: 



1, 



,6, 



(A5) 



where 77, ^, 7^ are Majorana fermions. This Majorana 
fermion representation enlarges the Hilbert space of the 
spins and generate a Z2 gauge redundancy on each siteQ- 
With this fermion representation, the Hamiltonian in 
Eq. (|A2p is given by 

3 



•K^exp[z((/)^ + (/>^0](A6) 

{pp') 



where m j. 



2a- 



n 



and exp(z(/)p) 



{jk)ep ^^jk 



(here we implicitly assumed that the order 



of {jk) in the product is taken counterclockwise). It is 
straightforward to show that the operators ^ii,i+e« are 
conserved in the enlarged Hilbert space, which makes the 
model exactly solvable since the Hamiltonian in Eq. ()A6p 
consists of quadratic fermions coupled with this static 
background Z2 gauge field Ui^iJ^^^ = ±1. The flux on 
each triangular plaquette (j)p is ±7r/2. 

For > and K » |J«|, |J4|, |J«|, |J4|, and |Jt|, 
the sum of neighboring (j)p and (j)p' must be zero in the 
ground state, which implies that the flux of up-triangle 
plaquettes must be opposite to the down-triangles. This 
flux pattern spontaneously breaks time reversal symme- 
try, as desired. So the Fermi surface, if exists, will be 
stable against weak perturbations, as shown below. In 
the following, we shall focus on the flux sector with 7r/2 
on up-triangles and — 7r/2 on down- triangles. 

We combine two Majorana fermions 77 and ^ on each 
site to form one complex fermions / with t] = f -\- and 

= ^(/^ ~ /)• Then the Hamiltonian Eq. (|A6p is given 

by 

+ [{J. + J'a) + i{Ja-J'J]fIfI+. 



h.c. 



(A7) 



The above Hamiltonian has an emergent U{1) symme- 
try (global Fermion number conservation) when Jq, = 
and Ja = —J'ai foi" which the pairing terms are identically 
zero. It is straightforward to obtain the single-particle 
dispersion of the hopping part: 

3 

ek = 4 ^ 1^ Jq, sin (k • Gq,) + Jq, cos (k • Gc^) —2Jj. (AS) 
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For = {0.3,0.3,0.3}, = {0.3,0.3,0.3}, = 
{1.0,1.0,1.0}, = {-0.8,-1.0,-1.2}, and Jj = 0.5, 
the Fermi surface of = is plotted in Fig. 2(b) It 
is clear that e_k 7^ for generic momentum k which 
is a consequence of broken time reversal and inversion 
symmetries. The feature makes the Fermi surface per- 
turbatively stable in the sense that adding weak pairing 
term cannot fully gap the Fermi surface. 

By Fourier transform, the full quadratic Hamiltonian 
is given by 



(A9) 



where ^l = {flf- 



<:) and 



At 



Ak 

-e-k 



(AlO) 



Here Ak = ELi - J^) + i{Ja + J^)] sm(k • e^.). 

The quasi-particle spectrum is given by 



^±(k) 



^^^-^ ± + |AkP. (All) 



The condition for E±{\i) = is given by 



(A12) 



When the system breaks both time reversal and inversion 
symmetries, e_k ^ Ck in general. Then, Eq. ()A12p is 
equivalent to only one real constraint. Since there are 
two momenta kx and ky, it is then clear that a line of 
zero-energy momentum points is generally allowed and 
stable. The Fermi surface of £'-(k) = is shown in 
Fig. 2(c)[ (If the system were time reversal invariant 
which requires ek = e_k and Ak is real; it is then clear 
that Eq. (|A12p would imply two constraints: ek = and 
Ak = 0, which generally allows at most nodal points, if 
not fully gapped.) 

A Z2 spin liquid with a stable pseudo-Fermi surface 
and broken time reversal symmetry is expected to exhibit 
finite thermal Hall conductivity which is linear in tem- 
perature. The thermal Hall conductivity is given bv(39| 



2 L.2 



klT 



3 h 



^ j d^kfxyO<.)ni^ 



(A13) 



where fxyO^) is the Berry curvature of the valence band 
in the momentum space and nk is one (zero) for occu- 
pied (empty) states. We denote /ik = E^=o^m(^)^^' 
where cr^ and a^'^'^ are identity and Pauli matrices, re- 
spectively. Then, the Berry curvature is given by 



Ixy 



(k) = — d(k) . dkJik) X dkjik) (A14) 



where d{k) = d{k)/\d{k)\ and d{k) 
{di{k)^d2{k)^ds{k)). It is clear that the Berry cur- 
vature cannot be zero everywhere when the phase 
of Ak is not uniform in the Brillouin zone. This is 
fulfilled when {Ja — J'^)/{Ja + Ja) ^^t all equal to 
one another. Ja = {0.3,0.3,0.3}, = {0.3,0.3,0.3}, 
Ja = {1.0,1.0,1.0}, = {-0.8,-1.0,-1.2}, and 
J7 = 0.5, we obtain the thermal Hall conductivity 

fc,,/r«o.2(^^). 



B. Estimate of pairing scale for nodal spin liquids 

For a Dirac node, the universal thermal conductivity 
is expected to be 



(Bl) 



where vp and v/\ are the velocities across and along the 
Fermi surface, respectively, at the nodes. Since 



3M 



0.4 X 



(B2) 



and the measured value is limT^o t^/T = 0.21^^ ^, 
we see that the velocity anistropy must be 



vf/va ^ 500, 



(B3) 



which is extremely large. 

For a pairing scaling A, we expect the velocity 
anistropy to be 



vf/va « J/ A ^ 500, 



(B4) 



which requires A ^ 0.5 K. Such an extremely small pair- 
ing scale relative to the exchange coupling J requires 
some explanation. It has been proposed that the pair- 
ing scale is roughly A ~ li^, which corresponds to a 
drop in the NMR relaxation rate. The small pairing scale 
has been suggested to be a conseqence of proximity to a 
nearby quantum critical point separating the nodal Z2 
spin liquid and the U{1) spin liquid. Above this scale, 
the system may be well-described by a U{1) spin liquid 
with spinon Fermi surface. A possible advantage of the 
pseudo Fermi surface proposals is that a small pairing 
scale is not required to explain the large thermal con- 
ductivity. Furthermore, if the ground state is indeed a 
nodal spin liquid, another natural explanation for the 
small pairing scale is a marginal instability of the Z4 
spin hquid with pseudo Fermi surface. 



